files = ['HH';'HV';'VH';'VV'];
pol = ['RL';'RR';'LL';'LR'];
Ts = [11,20,30,40,50,55,60,65,70,80,100,130,160,190,220,250,275,296];
figure
hold on
dataall=zeros(19,2958);
for num=1:length(Ts)
    filename=['VH',num2str(Ts(num)),'.txt'];
    data=readmatrix(filename);
    plot(data(:,1)/8.06554,data(:,2)+50*(num-1),'.-','Color',[(num)/20 0 0]);
    Es=data(:,1)/8.06554;
    dataall(num,:)=data(:,2);
    ylim([20,1000])
    title('VH')
end
figure
imagesc(Es,Ts,dataall)
colormap('jet')
clim([20,50])
figure
hold on
dataall=zeros(19,2958);
for num=1:length(Ts)
    filename=['HH',num2str(Ts(num)),'.txt'];
    data=readmatrix(filename);
    plot(data(:,1)/8.06554,data(:,2)+50*(num-1),'.-','Color',[(num)/20 0 0]);
    Es=data(:,1)/8.06554;
    dataall(num,:)=data(:,2);
    ylim([20,1000])
    title('HH')
end
figure
imagesc(Es,Ts,dataall)
colormap('jet')
clim([20,80])

figure


[X,Y]=meshgrid(Es,Ts)
mesh(X,Y,dataall)
colormap('jet')
clim([15,60])
xlim([55,70])